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Precise gap structure in LiFeAs (T c = 18K) given by ARPES studies offers us significant infor¬ 
mation to understand the pairing mechanism in iron-based superconductors. The most remarkable 
characteristics in LiFeAs gap structure would be that “the largest gap emerges on the tiny hole- 
pockets around Z point”. This result had been naturally explained in terms of the orbital-fluctuation 
scenario (T. Saito et ai, Pliys. Rev. B 90, 035104 (2014)), whereas an opposite result is obtained 
by the spin-fluctuation scenario. In this paper, we study the gap structure in LiFeAs by taking the 
spin-orbit interaction (SOI) into account, motivated by the recent ARPES studies that revealed the 
significant SOI-induced modification of the Fermi surface topology. For this purpose, we construct 
the two possible tight-binding models with finite SOI by referring the bandstructures given by dif¬ 
ferent ARPES groups. In addition, we extend the gap equation for multiorbital systems with finite 
SOI, and calculate the gap functions by applying the orbital-spin fluctuation theory. On the basis of 
both SOI-induced band structures, main characteristics of the gap structure in LiFeAs are naturally 
reproduced only in the presence of strong inter-orbital interactions between {d xz / yz -d xy ) orbitals. 
Thus, the experimental gap structure in LiFeAs is a strong evidence for the orbital-fluctuation 
pairing mechanism. 

PACS numbers: 74.70.Xa, 74.20.-z, 74.20.Rp 


I. INTRODUCTION 

The pairing mechanism of iron-based superconductors 
had been studied intensively as a central issue. Up to 
now, the spin-fluctuation-mediated s±-wave state with 
sign reversal [1-5] and orbital-fluctuation-mediated s++- 
wave state without sign reversal [6-8] had been proposed. 
According to theoretical studies for realistic multiorbital 
models, s±-wave state is fragile against impurities: s±- 
wave state with T c = 30K disappears when the residual 
resistivity po is as large as ~ 15 /rflcm for m*/mb ~ 3 
by 3d-transition metal impurities [9, 10]. Experimentally, 
the superconductivity in LaFeAsO+F, (K,Ba)Fe 2 As 2 and 
Ba(Fe,Co) 2 As 2 survives till p 0 ~ 300 /jQcm [11 -14], indi¬ 
cate that the s++-wave state occurs in these compounds. 
On the other hand, T c = 18K disappears only when 
po ~ 40 pUlcm in Ba(Fe,Ru) 2 As 2 by electron irradia¬ 
tion measurement [15]. In the neutron spectroscopy, the 
resonance-like peaks in various compounds are also ex¬ 
plained in terms of the s++-wave state by taking the re¬ 
alistic quasiparticle inelastic scattering into account [16]. 

In many Fe-based superconductors, the gap structure 
has been determined in detail by the angle-resolved pho¬ 
toemission spectroscopy (ARPES) studies. Although 
ARPES is not a sign-sensitive experiment, the detailed 
gap structure given by the ARPES offers us very use¬ 
ful information to distinguish the pairing mechanism. In 
fact, orbital dependence of the gap function is very cru¬ 
cial since it reflects the orbital dependence of the pairing 
interaction strength. In the spin-fluctuation mechanism, 
the superconductivity is mainly induced by the intra¬ 
orbital scattering of the Cooper-pairs, since the spin- 


fluctuations are driven by the intra-orbital Coulomb in¬ 
teraction U. For this reason, the gap magnitude |A(fe)| 
tends to largely depend on the orbital character of the 
Fermi surface in the spin-fluctuation mechanism. In con¬ 
trast, in the orbital-fluctuation mechanism, the super¬ 
conductivity is induced by the inter-orbital scattering of 
the Cooper-pairs. For this reason, the orbital dependence 
of |A(fe)| tends to be moderate in the orbital-fluctuation 
mechanism. 

Motivated by the above discussion, the gap magni¬ 
tude on the d z 2 -orbital hole-pocket around Z point in 
Bal22 and Srl22 compounds had been studied inten¬ 
sively. Because of the absence of the intra-d z 2 -orbital 
nesting, the spin-fluctuation scenario predicts the emer¬ 
gence of the ’’horizontal node” on the hole-pocket around 
Z point [17, 18]. In contrast, the horizontal node is 
absent in the orbital-fluctuation-mediated s-wave state, 
since the d z 2 -orbital electrons contribute to the orbital- 
fluctuations [18]. In BaFe 2 (As,P )2 and Ba(Fe,Co) 2 As 2 , 
several ARPES groups reported the absence of the hor¬ 
izontal node [19-21], whereas horizontal node was re¬ 
ported in Ref. [22]. 

For studying the pairing mechanism, LiFeAs is an suit¬ 
able compound since the detailed band structure and the 
superconducting gap functions had been determined by 
several ARPES groups [23, 24]. LiFeAs is an ideal sys¬ 
tem for ARPES studies in that the charge neutral cleav¬ 
age plane exists and clean single crystals can be synthe¬ 
sized. Therefore, the intrinsic gap structure free from 
extrinsic effects (such as the impurity effect) is obtained 
by ARPES studies. The experimental band structure 
and Fermi surfaces (FSs) are shown in Fig. 1: Reflect- 
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ing the bad nesting between the hole-like FSs (h-FSs) 
and electron-like FSs (e-FSs) in LiFeAs, moderate spin- 
fluctuations are observed by NMR studies [25] and neu¬ 
tron scattering studies [26-28]. Interestingly, in LiFeAs, 
the largest gap appears in the tiny three-dimensional 
(3D) hole-pockets (h-FSl,2) according to the ARPES 
studies [23, 241 and the Scanning Tunneling Microscopy 
(STM) study [29]. 




FIG. 1: (Color online) (a) The dispersion of the band struc¬ 
ture between points A and B for A = 0.05 eV (red lines) 
and 0 (blue lines), respectively. We call the former the “SOI- 
induced bandstructure (I)”. A and B points are shown in (b). 
The corresponding FSs for (b) A = 0 and (c) 0.05 eV in the 
k z = 7r plane are shown respectively. In (c), h-FSl disappears 
due to the SOI. Green, red, and blue lines correspond to xz, 
yz, and xy orbitals, respectively. 

Based on the tight-binding model for LiFeAs con¬ 
structed from the ARPES band structure, some au¬ 
thors studied the possible pairing scenarios. In the spin- 
fluctuation mechanism for the s±-wave state, the gaps 
on the tiny hole-pockets, h-FSl,2 in Fig. 1 (b), are very 
small, because of the bad intra-orbital nesting between 
h-FSl,2 and e-FSs [31, 32]. In contrast, in the orbital- 
fluctuation mechanism for the s_|_+-wave state, the gaps 
on the tiny hole-pockets are the largest, because of the 
good inter-orbital nesting [32]. The latter theoretical re¬ 
sult is consistent with experiments [23, 24, 29]. 

However, although the SOI had been neglected in pre¬ 
vious theoretical studies in Refs. [31, 32], recent ARPES 
studies for LiFeAs revealed that the SOI is never negligi¬ 
ble to construct the FSs and the band structure near 


the Fermi level. For instance, the degeneracy of the 
d X z/ yz hole bands on the F-Z line is lifted by the SOI 
[33-35]. For this reason, one of the tiny two hole-pockets 
(h-FSl,2) disappears due to the SOI, as shown in Fig. 1 
(c) (or Fig. 6 (d)). Similar disappearance of h-FSl due 
to SOI had been observed in FeSe [36-41]. Since the tiny 
hole-pockets play important roles in the pairing mecha¬ 
nism, a theoretical study for LiFeAs including the SOI 
had been highly required. 

In this paper, we study the gap structure of LiFeAs by 
taking the effect of the SOI into account, in order to dis¬ 
tinguish the correct pairing scenario. For this purpose, 
we construct two possible 3D ten-orbital models, which 
correspond to the different experimental band structures 
reported in Refs. [33, 34]. Both models possess a sin¬ 
gle tiny hole-pocket around Z point. We also introduce 
both the multiorbital Coulomb interaction (U, U', J) and 
the quadrupole interaction ( g ): The latter will origi¬ 
nate from the Aslamazov-Larkin type vertex correction 
(AL-VC), which describes the many-body effect beyond 
the mean-field approximation [7, 8, 42, 43]. By us¬ 
ing the random-phase-approximation (RPA), the orbital- 
and spin-fluctuations develop with increasing g and U, re¬ 
spectively. Then, we extend the gap equation for multior- 
bital systems with finite SOI, and calculate the gap func¬ 
tion by using the pairing interaction given by the RPA. 
In the s++-wave state due to orbital-fluctuations, the gap 
on the smallest h-FS is the largest of all. In contrast, in 
the s±-wave state due to spin-fluctuations, the gap on 
the smallest h-FS is almost zero. Therefore, the super¬ 
conductivity in LiFeAs originates from the inter-orbital 
idxz/yz-dxy) interactions due to orbital-fluctuations, not 
from the intra-orbital spin-fluctuation interactions. 

When both the orbital- and spin-fluctuations develop 
comparably, we frequently obtain the “hole-s±-state”, in 
which the gap structure has the “sign reversal between h- 
FSs”. This interesting gap structure had been discussed 
in (Ba,K)Fe 2 As 2 experimentally [44, 45]. In Appendix 
A, we show that the hole-s±-wave gap structure can be 
obtained in the 3D model for Bal22 system. 


II. FORMALISM 

A. SOI-induced bandstructure (I) and interaction 
terms 

In this paper, we set x and y axes parallel to the nearest 
Fe-Fe bonds, and the orbital z 2 , xz, yz, xy, and x 2 — y 2 
of Fe-A (B) site are denoted as 1, 2, 3, 4, and 5 (6, 7, 
8, 9, and 10) respectively. We used the 3D ten-orbital 
tight-binding model which is given in Ref. [31]. This 
model is obtained by fitting the experimentally observed 
dispersion reported in Ref. [23]. Note that the band 
renormalization due to the mass enhancement m*/mb = 
2 ~ 3 is taken into account in the present model. 

We stress that the experimental FSs of LiFeAs in Fig. 
1, on which the present study is based, are very differ- 
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ent from the FSs given by the density functional theory [46] 
(DFT), in which li-FSl,2 predicted by the DFT are much 
larger. Better agreement between theory and ARPES is 
achieved by the LDA+DMFT study [30], since the FS1,2 
shrinks due to the orbital-dependent self-energy that is 
absent in the LDA. 

In the absence of the SOI, the kinetic term of the ten- 
orbital model is given as , 


X S (q) = 


X° (g) 

i - f s x° {q) ’ 


X C (Q) 


X° (g) 

i-f c (wi)x° (q)' 


( 4 ) 

( 5 ) 


H° = ^2 tl,m{ R a ~ R b)cl a (Ra)CmA R b) 

ablmcr 

= E \ T, t ^( R a)e ikA Acl(k)c mcr (k) (1) 

where l,m = 1 — 5 (6 — 10) represent the d orbitals at 
Fe-A (Fe-B) atom, and a = ±1 is the spin index. R a is 
the Fe-site position measured from an Fe-A atom, c' la .(k) 
is the creation operator of the d electron, and ti <m ( R ) is 
the hopping integral. The values of ti >m ( R ) are shown in 
Ref. [31]. 

Next, we introduce the interaction terms. We con¬ 
sider both the multiorbital Coulomb interaction (£7, £/', 
J = (£7 — U')/ 2) and the quadrupole interaction. The 
quadrupole interaction Hamiltonian is given as 


(r S )ql,i 3 Z4 — 'W x { 


f U, l\ — I 2 — I 3 — h 

U', h = I3 ^ h = U 

T, l\ — I 2 / ^3 — ^4 

J', h = I 4 ^ h = h 

[0, otherwise 


( 6 ) 


f c ( Wi ) = -C-2V quad {oj l ), 


( 7 ) 


(C)qf 2 ,j 3 Z4 — 'W x ^ 


(£7, 1 1 — I 2 — h — h 

— U' + 2 J, l\ = ?3 ^ I 2 = I 4 

2 U' — J, h = h ^ h = I 4 (8) 

J', h = I 4 ^ l 2 = I 3 

0. otherwise 


site 

Kjuad = E (o a yz ■ d a yz + o a xz ■ O a xz ) (2) 

a 

where g(wi) = gu^/^f + Wq), and g = g(0) is the 
quadrupole interaction at u>i = 0. u >0 is the cutoff en¬ 
ergy of the quadrupole interaction. Op = [T/ Op m“ m 

(m“ m = Xlcr c L(- R a) c ^(- R “)) is the quadrupole oper¬ 
ator at site R a introduced in Ref. [6]: The non-zero 
coefficients of Op" 1 = o™’* are o x ’ z = o x ’ z = V?>o l x ' z — 1, 
and — Oy’ z = Oy’ z = \/3o y ’ z = 1. Thus, V qua d has many 
non-zero inter-orbital elements. As explained in Ref. 
[6], g is given by in-plane Fe-ion oscillations. Also, the 
AL-VC due to Coulomb interaction gives large effective 
quadrupole interaction g [7]. 

Now, we analyze the present model by applying the 
RPA. The 32 x 32 x 16 k meshes are used in the numerical 
study. We fix the temperature at T = 0.01 eV, the cutoff 
energy at ujq = 0.02 eV, and set the filling of each Fe-site 
as n = 6.0. The irreducible susceptibility is given as 

Xw,mm' ( 9 ) = ~ jy E Gl,m (& + d) Gm',V (&) j. (3) 

k 

where q = (q,u>i) and k = (k, e n ). e n = (2 n + l)nT 
and uii = 2lnT are the fermion and boson Matsubara 
frequencies. G(k) = [ie n + g — h°(fc)] _1 is the Green 
function in the orbital representation, where h°(k) is the 
matrix representation of H° for the momentum k, and 
g is the chemical potential. In the RPA, the spin and 
charge susceptibilities in the matrix form are given by 


a, /3(= A, B ) is Fe-site. In this paper, we set U = £7'+2 J, 
J = J', and J = £7/6. 

In the RPA, the enhancement of the spin susceptibility 
X s originates from the intra-orbital Coulomb interaction 
U and the “intra-orbital nesting” of the FSs. On the 
other hand, the enhancement of \ c originates from the 
quadrupole-quadrupole interaction in Eq. (2) and the 
“inter-orbital nesting” of the FSs. The magnetic (orbital) 
order occurs when ct s ( c ) = 1, where a s ( c ) is the the maxi¬ 
mum eigenvalue of T s ^ c - , x^°' ) (<?, 0), called the spin (charge) 
Stoner factor. Here, the magnetic order is realized when 
U reaches E7 cr = 0.448 eV. Also, the orbital order is re¬ 
alized when g reaches g c r = 0.132 eV for £7 = 0 [32], 
The used interactions are smaller since the energy-scale 
had been renormalized by z = (m*/mb)^ 1 in the present 
tight-binding model. 

B. Spin-orbit interaction 

In this subsection, we introduce the SOI. In the pres¬ 
ence of SOI for d electron, the 20 x 20 matrix expression 
of the total Hamiltonian at k is given by 

h(k) = fh°(k) + \l z /2 \{l x - il y )/2 \ 

1 ; \\(i x + ii y )/2 h 0 (k)-\i z /2j 

htt(fc) §t4(fe)\ 
h*A(k) /rA(fc )J 

where the first and the second rows (columns) correspond 
to t-spin and 4_-spin. A is coupling constant of SOI. The 
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matrix elements for l for d-orbital, which is diagonal with 
respect to site, are given by [47] 


( 0 

0 

i 0 

1 0N \ 

0 

0 0 

i 

0 


O 

O 

0 

i —i 

0 

-i 0 

0 

i 0 

V o 

0 i 

0 

i o) 

j 0 

-V3i 0 

0 

°\ 

V3i 

0 0 

0 

—i 

0 

0 0 

—i 

; o 

0 

0 i 

0 

0 

V o 

i 0 

0 

°/ 

/0 0 

0 0 



0 0 

—i 0 

0 


0 i 

0 0 

0 


0 0 

0 0 

2 i 


\0 0 

0 -2 i 

oj 



( 10 ) 


( 11 ) 


( 12 ) 


where the first to fifth rows (columns) correspond to d 
orbitals z 2 , xz, yz, xy, x 2 — y 2 , respectively. 

The red (blue) lines in Fig. 1 (a) shows the band struc¬ 
ture of LiFeAs for A = 0.05 eV (A = 0), respectively. 
Points A and B are shown in Fig. 1 (b). When A = 0, 
the two bands composed of d xz + d yz orbital near the 
Fermi level at Z point are degenerated, however, this de¬ 
generacy is lifted by the SOI, which seems to be observed 
by recent ARPES measurements in LiFeAs [33, 34] as 
well as in FeSe [36-41]. Figures 1 (b) and (c) show the 
FSs of LiFeAs for A = 0 and 0.05 eV, respectively. For 
this reason, h-FSl disappears by introducing the SOL 
This model is consistent with the ARPES measurement 
by Miao et al. in the way that only outer band of two 
d xz /y Z -orbital hole band at Z point crosses the Fermi level 
[34]. In Sec. IV, we also construct the tight-binding 
model that gives the ARPES band structure reported by 
Borisenko et al [33]. 

Since the energy-scale of the present tight-binding 
model had been renormalized by the factor z = 
(?n*/TOfc) _1 , we should use a renormalized A in Eq. (9) 
in principle. However, we hereafter use a typical A for Fe 
ion (A = 0.05 eV) [47] without considering the renormal¬ 
ization, in order to clarify the effect of the SOL In fact, 
we had verified that the obtained gap structures in later 
sections are unchanged by using smaller SOI coupling, 
say A ~ 0.02 eV. 


C. Gap equation 


given as 


Ae A ss (fc,e„) — 


7FT 

(2ir) 3 


FS 


EEE 


M A 


f dk' FSM 

FS M V M (k') 

Af A (fc',e m ) 


xViH- A (k,k',e n -em) AA ,r , (13) 


where Ae is the eigenvalue of the gap equation. Ae 
reaches unity at T = T c . We solve the gap function for 
the largest eigenvalue, since it corresponds to the pairing 
state for the highest T c . L and M denote the FSs, and 
E and A denote the pseudo-spin (f|-, JJ.) in the presence of 
the SOL Note that E = — E. A£g(fe,e n ) is the pseudo¬ 
spin singlet gap function on the L-th FS. In Eq. (13), we 
perform the surface integral on the M -th FS. The paring 
interaction V in Eq. (13) for the pseudo-spin singlet pair 
is given by 


E£AA( fc ’ fe, ’ e « “ e m) 

= E E U£*(k)Uffi*(-k) 

Imm'l' crAA'cr' 

x Vimm'v (k — k , e n — e m ) 

xU^(-k')U^Uk'), (14) 


where a, A mean real spin (4_,t) and t/^ s (fc) = 
(fe; la\k; HE) is the unitary matrix connecting between 
the band representation and the orbital one. Here, the 
pairing interaction in the orbital basis is given by 


rcrXX <j 


Vlmm'l'dcr + Vlmm'l l<T cr\ ' & V \' 

| ^ Imm'l ' ^ Imm'l'i 

b 

II 

II 

-c 

II 

b 

i ^ Imm'l ' ^ Imm'l' 

<r = A ^ A' = a 1 

1 ‘ZVlmm'l’i 

u = A' A = er' 


otherwise, 


(15) 


v c = E c x c f c , V s = if s x s f s . (16) 


In this paper, V? (£ = Cj S J i s calculated without the SOI, 
since it is verified that the SOI is negligible except at the 
vicinity of the magnetic critical point (a s « 1). 

Now, we introduce the relation between f7(^ s (fc) and 
UiL{—k) to preserve the time-reversal symmetry be¬ 
tween the states |fc;LE) and | — fe;LE). The 20 x 20 
Hamiltonian h(k) is given in Eq. (9). Then, h(—k)* is 
given by 


In this subsection, we introduce the linearized gap 
equation in the presence of the SOL In the 3D model, 
it is difficult to analyze fine momentum dependence of 
the SC gap on the FSs by using a conventional 3D fe- 
meshes. In this paper, we used 48 x 16 fe points for each 
FS sheet in solving the gap equation in order to obtain 
the precise gap structure [48]. In the absence of impu¬ 
rities, the linearized gap equation for the singlet state is 


/feo(-fc)* + Af*/2 \(i x -ii y )*/2 \ 

V A(l + z4)*/2 feo(-fe)* - \i* z /2j 

(h 0 (k)-\l z /2 —\{l x -\-iiy)/2\ 

\ — A (l x ~ ily)/2 ho(k) + Xl z /2 ) 
f feH(fc) —fell(fe)\ 

\-h^(k) fctt(fc) ) ‘ 
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Here, we express the unitary matrix as U(k) 
^l(fcj) that diagonalize h(k). That is, 


III. SUPERCONDUCTING GAP FOR THE 
SOI-INDUCED BANDSTRUCTURE (I) 


U(k)'h(k)U(k) = ) , (18) 

where E(k) is diagonal. By using Eqs. (17) and (18), we 
obtain that 

( UW(k)i -U™{k) t\ u 
\~u^{k )t cm(fc)t ) [ ] 

U i m ~U u (k)\ (E{k) 0 \ 

-U^{k) U™{k) ) V6 E(k)J ' 


In this section, we numerically analyze the linearized 
Eliashberg equation, Eq. (23), on the basis of the SOI- 
induced bandstructure in Fig. 1 (a), which we call the 
“bandstructure (I)”. We show the obtained 3D gap func¬ 
tion A L (d,k z ) on the FS sheet L. Here, we divide the 
variables 9 = [0,27r] and k z = [— n, 7r] into 48 and 16 
meshes, respectively, and use 512 Matsubara frequencies. 
The pairing interaction in Eq. (14) is given by the RPA. 


A. Orbital-fluctuation-mediated s++-wave state 


Therefore, U(—k) is related to U(k) as 


U{-k) 


( U^(k)* -U^{k)*\ 
\-u^(k)* cr^(fe)* )' 


( 20 ) 


That is, the relations between Uf^{k) and Uf^{—k) are 
summarized as 

u!H-k) = uH(kr, 

U${-k) = u™(k)*, 

u!h-k) = -utfik)*, 

U^{-k) = -U^{k)*. ( 21 ) 

These relations are necessary to obtain physical gap func¬ 
tions. 

Next, we transform the gap equation (13) to a more 
convenient expression. By setting E =-f( (E =JJ.) and 
taking the summation of A, gap equation is given as 


A E A^(fc, e n ) (27r)3 E E J FSM v M {kf) 


FS 


rlk' 

afi, FSM 


X \ k ', €n ~ e rn) 


A ¥Jk',e m ) 


+ k !, e n — e m ) 


A fJk',e m )' 


( 22 ) 


We first discuss the s ++ -wave state that is realized 
by orbital-fluctuations in the presence of the SOI: Fig¬ 
ure 2 shows the obtained gap functions in the case of 
g = 0.129 eV and U = 0 (a c = 0.98) in the k z = tt 
plane. The overall gap structure is essentially unchanged 
for a c = 0.90 ~ 0.98 [32]. As for the hole pockets, the 
gap functions on the h-FS2 composed of ( d xz ,d yz ) or¬ 
bitals are the largest in magnitude, while the gap on the 
h-FS3 composed of the d xy orbital is the smallest. This 
result is quantitatively consistent with the experimental 
reports in Refs. [23, 34]: Experimental gap functions in 
Ref. [23] are shown by dots in Fig. 2. This result is also 
similar to our previous result without SOI, regardless of 
the presence or absence of the h-FSl [32]. Here, we ad¬ 
just the magnitude of gap functions given by solving the 
linearized gap equation. 


h-FS2 ( xz/yz ) 


h-FS3 (xy) 



Since the relation A^(fc) = — A^(/c) is satisfied, the 
gap equation for the singlet pairing state is given as 


A e A L {k,e n ) = 


7 tT 


FS 


EE 


rlk’ 

U *‘FSM 

M(UI\ 


M 


FSM V AI (k') 


x { V imv(k,k',e n - e m ) 

(23) 


where A L {k,e n ) = A^(fe,e„). 


FIG. 2: (Color online) Obtained gap functions for [7 = 0 and 
g = 0.129 eV (s++-wave state) in the k z = 7r-plane for the 
SOI-induced bandstructure (I). The dots represent the experi¬ 
mental data (in unit meV) given by the ARPES measurement 
in Ref. [23]. 

As for the electron pockets, the gap of the inner e-FS 
is larger than that of the outer e-FS, and has the local 
maxima at 9 = 0 and n/2 and the minima at 9 = t r/4. 
This result is also consistent with the experimental data 
[23], and it is very similar to our previous result without 
SOI [32]. The obtained gap structure is almost indepen- 
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dent of k z except for the three dimensional pocket h-FS2. 
These obtained results are essentially independent of the 
strength of g. 

Next, we discuss the origin of the orbital and FS de¬ 
pendencies of the gap functions. For finite g , orbital- 
fluctuations develop due to the good inter orbital nesting 
between h-FS2 (orbital 2,3) and e-FSs (orbital 4). For 
this reason, the maximum gap is realized on h-FS2 (All), 
and inner e-FS (A® ut ) at 9 = 0 and 7 t/ 2. This result is not 
changed by considering the SOL Therefore, the experi¬ 
mentally observed gap functions are understood based on 
the orbital-fluctuation theory. 


B. Spin-fluctuation-mediated s±-wave state 

Next, we discuss the s±-wave state realized by spin- 
fluctuations: Figure 3 shows the obtained gap structure 
in the case of g = 0 and U = 0.439 eV (a s = 0.98) in 
the k z = 7r plane. The overall gap structure is essentially 
unchanged for a s = 0.90 ~ 0.98 [32]. The gap function 
on each FS is almost independent of k z . The obtained 
gap structure is essentially independent of the value of 
U. The obtained very small gap on the h-FS2 is totally 
opposite to the experimental data shown by dots. 


5 

<1 

0 

o n/2 0 nil 

o 9 

FIG. 3: (Color online) Obtained gap functions for U = 0.439 
eV and g = 0 (s±-wave state) in the k z = 7r-plane. The dots 
represent the experimental data (in unit meV) [23]. 



C. Coexistence of orbital- and spin-fluctuations 

Now, we discuss the superconducting state when the 
orbital- and spin-fluctuations coexist. In the previous 
results without SOI, coexistence of orbital- and spin- 
fluctuations leads to an exotic s-wave state, which has 
sign reversal within h-FSs: A|Ag < 0 [32]. Here, we 
show that this ”hole-s±-state” is also realized in the pres¬ 
ence of the SOI. 

Figure 4 (a) shows the obtained gap functions in the 
case of g = 0.125 eV and U = 0.200 eV. The realized 
Stoner factors are a c = 0.98 and a s = 0.45. In this 
case, the orbital-fluctuations are much larger than the 
spin-fluctuations, and therefore we obtain the s + _|_-wave 
state. Except for h-FS3, the obtained gap structures are 
similar to those of the ’’pure s ++ state” without U in 
Fig. 2. Due to the moderate spin-fluctuations on the d xy 
orbital, the anisotropy of A| is enlarged. 



In addition, the obtained 9 dependence of the gap on 
the e-FSs and h-FS3 is very different from the experimen¬ 
tal data. Both the gap on the h-FS3 and e-FSs show the 
maximum values at 9 = 7t/ 4, since they are connected 
by the wave vector of spin-fluctuations Q ~ (n, 0), (0,7r). 
Moreover, the gap function on h-FS3 has eight nodes, 
which is inconsistent with experiments. These results 
are very similar to our previous results without SOI [32]. 
Note that these eight nodes on the hole-FSs disappear in 
the case of J/U = 0.4 (U' = U — 2 J = 0.217) in accor¬ 
dance with the previous study in Ref. [31]. 


FIG. 4: (Color online) (a) Obtained gap functions for U = 
0.200 eV and g = 0.125 (s++-wave state) in the k z = 7r-plane. 
The dots represent the experimental data (in unit meV) [23]. 
(b) Obtained gap functions for U = 0.380 eV and g = 0.122 
(hole-s±-wave state) in the k z = 7r-plane. 

If we increase the value of U further, we obtain the 
hole-s±-state with the sign reversal within hole pockets. 
Figure 4 (b) shows the obtained gap functions in the 
case of g = 0.122 eV and U = 0.380 eV (a c = 0.98 and 
a s = 0.85). Here, only A| is negative. In this hole-s±- 
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wave state, the gap structure of Ajj and that of the e-FSs 
(A e ) are qualitatively similar to those in the s ++ -wave 
state in Fig. 2. On the other hand, A| becomes very 
anisotropic, similarly to A| in the s±-wave state in Fig. 
3. 

In the hole-s±-wave state given by the coexistence of 
orbital- and spin-fluctuations, Ajj • A e is positive due to 
orbital-fluctuations, whereas A 3 • A® is negative due to 
spin-fluctuations. Although the gap on the tiny hole- 
pocket in the hole-s±-wave state takes the largest value, 
the overall gap structure is not consistent with ARPES 
measurement in Ref. [23], The present hole-s±-wave 
due to cooperation of orbital- and spin-fluctuation would 
be realized in other iron based superconductors. For 
instance, the hole-s±-wave state is first discussed in 
Bai_ x K a ;Fe 2 As 2 based on the thermal conductivity and 
penetration depth measurements [44], in addition to the 
recent ARPES study [45]. 



FIG. 5: (Color online) a c -a a phase diagram of the gap struc¬ 
ture in LiFeAs obtained for the bandstructure (I) in Sec. III. 
Similar phase diagram is obtained for the bandstructure (II). 
The gap functions at each point a ~ c are shown. 

Figure 5 shows the obtained a c -a s phase diagram of 
the gap structure in LiFeAs for A = 0.05 eV (bandstruc¬ 
ture (I)). As expected, the s±-state (s++-state) is real¬ 
ized for a wide region of a c < a s (a c > a s ). When both 
a c and a s are close to unity, we obtain the hole-s± state 
gap in a wide region. The gap structure at each point a 
~ c is shown in the figure. In the region ” A| ~ 0”, the 
obtained A|(0) is very small and has line nodes, whereas 
the gaps on the other FSs are large and their signs are 
the same. Thus, nearly s ++ -wave state is realized. Thus, 
various types of s-wave gap structures are realized due to 
the coexistence of orbital- and spin-fluctuations. 

Finally, we comment that the hole-s± state is also re¬ 
alized in the presence of strong repulsive pairing inter¬ 
action within the h-FSs, as discussed in Refs. [49, 50]. 
The authors in Ref. [51] also discussed the hole-s± state 
due to strong repulsive interaction between d xy orbital 


in h-FS3 and d XZtVZ orbitals in e-FSs [51]. In the lat¬ 
ter scenario, the gaps on e-FSs are nodal in the presence 
of the SOI. In contrast, the hole-s± state in the present 
paper is realized by the cooperation between the “attrac¬ 
tive interaction among (d xz ,d yz )- and d X 2 / -orbitals” and 
“repulsive interaction on the d^y-orbital”. 


IV. ANOTHER SOI-INDUCED 
BANDSTRUCTURE AND GAP FUNCTION FOR 
LIFEAS 

A. SOI-induced bandstructure (II) 

Recently, the 3D bandstructure in LiFeAs had been 
carefully measured by ARPES study in Ref. [33]. The 
authors observed the strongly dispersing d z 2 Td^-orbital 
band crossing the Fermi level along T-Z line, with the 
minimum at Z point. Although this dispersion is ab¬ 
sent in bandstructure (I), it can be reproduced by low¬ 
ering the dcy-orbital level below the Fermi level at the 
Z point without changing d^-orbital h-FS3. To change 
the dj-orbital level by A£) (l = xz,yz,xy,x 2 — y 2 ,z 2 ) 
only around Z point, we introduce the additional intra¬ 
orbital hopping integrals shown in TABLE I. In this pro¬ 
cedure, the dj-level at M point in the unfolded Brillouin 
zone is unchanged. Figure 6 (a) shows the bandstructure 
for A E xy = —0.6 eV and A E xz = A E yz = 0.01 eV in 
the absence of the SOL Due to large negative A E xy , the 
strongly dispersing band along the F-Z line, shown by 
brown line, is reproduced. The energy of this dispersion 
is +0.21 eV (-0.51 eV) at T (Z) point. The Fermi level is 
set as zero. The FSs in the k y = 0 plane are shown in Fig. 
6 (b): After introducing large negative A E xy , d xy -leve\ 
is lower than d xz / yz -\eve\ at Z-point, and the former level 
(latter level) corresponds to the energy of the hole-band 
for h-FSl (h-FS2) at Z-point. As a result, the degen¬ 
eracy of the hole-bands for h-FSl,2 at Z point is lifted, 
and only h-FS2 remains in Figs. 1 (a)-(b) whereas the 
hole-band for h-FSl is below —0.1 eV. Thus, the number 
of the tiny 3D hole-pockets is reduced to one. 

In the next stage, we introduce the SOI (A = 0.05 
eV): The realized SOI-induced bandstructure is shown 
in Fig. 6 (c), which we call the “SOI-induced bandstruc¬ 
ture (II)”. Due to the SOI, the degenerated d xz / yz - orbital 
hole-bands along the T-Z line, is split into two bands, as 
shown by purple lines in Fig. 6 (c). As a result, the size 
of h-FS2 is reduced as explained in Ref. [33] in detail. 
(In contrast, the higher band forms the tiny hole-pocket 
in the SOI-induced bandstructure (I) in Sec. III.) 

The obtained SOI-induced FSs in the k y = 0 plane are 
shown in Fig. 6 (d): In addition to the tiny and large 
hole-pockets h-FS2 and h-FS3, there are fine h-FSs due 
to the Dirac dispersion at \k z \ ~ 0.5-7T on the fcj-line. The 
overall bandstructure near the Fermi level as well as the 
SOI-induced change in the topology of the FSs reported 
in Ref. [33] are qualitatively reproduced in the present 
bandstructure (II). 















TABLE I: Additional hopping integrals A ti,i(R) to change the dj-orbital level by AEi only around the Z point, without changing 
the dj-level at M point in the unfolded Brillouin zone, si = +1 ( — 1) for l = z 2 ,xy,x 2 — y 2 ( xz,yz). 


R 

[0,0,01 

[1,0,0] 

[1,1,0] [2,0,0] [2,1,0] 

A t,,i(R) 

(1/16) AEi 

—(3s;/64)A£) 

(1/32)A£) (1/64)A£) -(s;/128)A£) 

R 


[0,0,1] 

[L0,1] 

[1,1,1] [2,0,1] [2,1,1] 

A t u (R) 

- 

(1/32)A Ei 

(3si/128)A£) 

— (1/64)A E l —(1/128)A E t (s t /256)AE l 


B. Superconducting gap for the bandstructure (II) 

Next, we analyze the linearized Eliashberg equation, 
Eq. (23), for the SOI-induced bandstructure (II): The 
FSs in the k z = it plane are shown in Fig. 7 (a). We 
stress that the orbital character of h-FS2 in Fig. 7 (a) is 
essentially the same as that of h-FS2 in Fig. 1 (c). This 
fact indicates that similar gap structures are obtained 
based on both bandstructures (I) and (II). 

The pairing interaction is given by Eqs. (14) and (15) 
in the present model with A£). We first discuss the s++- 
wave state realized by orbital-fluctuations: Figure 7 (b) 
shows the obtained gap functions in the case of g = 0.127 
eV and U = 0 ( a c = 0.98) in the k z = it plane. As for the 
hole pockets, the gap functions on the li-FS2 composed 
of (d xz , d yz ) orbitals are the largest, while the gap on the 
h-FS3 composed of the d xy orbital is the smallest. These 
results are quantitatively consistent with the experimen¬ 
tal data [23] shown by dots. As for the electron pockets, 
the gap of inner e-FS is larger than that of outer e-FS, 
and their 9 dependences are qualitatively consistent with 
experimental data in Ref. [23]. The obtained overall 
gap structure is similar to that for the bandstructure (I) 
shown in Sec. Ill A, regardless of significantly different 
bandstructures. 

Next, we discuss the s±-wave state realized by spin- 
fluctuations: Figure 7 (c) shows the obtained gap struc¬ 
ture in the case of U = 0.434 eV and g — 0 (a s = 0.98) in 
the k z = 7T plane. The obtained very small gap functions 
on the h-FS2 is opposite to the experimental data [23] 
shown by dots. In addition, the obtained 6 dependences 
of the gaps on the e-FSs and h-FS3 are very different from 
the experimental data: Both the gaps on the h-FS3 and 
e-FSs show the maximum values at 9 = 7r/4 since they 
are connected by the wave vector of spin-fluctuations 
Q ~ (it, 0), (0,7r). These results are similar to those for 
the bandstructure (I) shown in Sec. IIIB. 


V. SUMMARY 

In this paper, we studied the pairing mechanism in 
LiFeAs by taking the SOI into account. For this purpose, 
we constructed possible two realistic 3D tight-binding 
models with SOI, according the experimental band struc¬ 
tures [33, 34]. In both SOI-induced bandstructures, one 
of the 3D tiny hole-pockets around the Z point disappears 
due to the SOI. The orbital character of the realized tiny 


single hole-pocket is same in both models. Experimen¬ 
tally, the superconducting gap takes the largest value on 
the tiny single hole-pocket. 

Next, we analyzed the 3D gap structure in order to 
determine the correct pairing mechanism. In the orbital- 
fluctuation scenario for the s ++ -wave state, the pairing 
interaction is given by the orbital-fluctuations driven by 
the inter-orbital nesting between ( d xz / yz , d xy ) orbitals. 
The obtained gap function on the tiny hole-pocket is the 
largest, by reflecting the better inter-orbital nesting be¬ 
tween the tiny hole-pocket and the electron-pockets. The 
obtained numerical results shown in Fig. 2 and Fig. 7 (b) 
are quantitatively consistent with the experimental gap 
structure in LiFeAs. In contrast, in the spin-fluctuation 
scenario for the s±-wave state, the obtained gap function 
on the tiny hole-pocket is the smallest, by reflecting the 
bad intra-orbital nesting between the tiny hole-pocket 
and the electron-pockets. Thus, the experimental gap 
structure in LiFeAs is unable to be explained by this sce¬ 
nario, as demonstrated in Fig. 3 and in Fig. 7 (c). 

Therefore, it is concluded that the inter-orbital 
(d xz / yz -d xy ) interactions, which are given by orbital- 
fluctuations, cause the superconductivity in LiFeAs. The 
present study reinforces the results obtained in our pre¬ 
vious study for LiFeAs without the SOI [32]. When 
both the orbital- and spin-fluctuations develop compa¬ 
rably, we frequently obtain the “hole-s±-state”, in which 
the gap structure has the “sign reversal between h-FSs”. 
This interesting gap structure had been discussed in 
Bai_ x K x Fe 2 As 2 based on the thermal conductivity and 
penetration depth measurement [44]. In Appendix A, 
we show that the hole-s±-wave gap structure can be ob¬ 
tained in the model of Bal22 system. 
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Appendix A: Hole-s±-wave state due to coexistence 
of orbital- and spin-fluctuations in Bal22 systems 

In this appendix, we show that the hole-s±-wave 
can be realized in Bal22 compounds under the coexis¬ 
tence of orbital- and spin-fluctuations. To construct the 
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(a) X=OeV 




k y 


FIG. 6: (Color online) (a) The band structure of the ten- 
orbital model for A E xy = —0.6 eV and A E xz = A E yz = 0.01 
eV without the SOI (A = 0) and (b) the FSs in the k y — 0 
plane for A = 0. Due to large negative A E xy , h-FSl around Z 
point disappears, (c) The “SOI-induced bandstructure (II)” 
with A = 0.05 eV, and (d) the corresponding FSs in the k y = 0 
plane. The size of h-FS2 is reduced by introducing the SOI. 


bandstructure, we perform a local-density-approximation 
(LDA) band calculation for BaFe 2 As 2 using the WIEN2K 
code based on the experimental crystal structure. Next, 
we derive the ten-orbital tight-binding model using the 
WANNIER90 code and the WIEN2WANNIER inter¬ 
face [52]. However, the d^-orbital h-FS given by the 


(a) 4= if plane A.=0.05eV 



FIG. 7: (Color online) (a) FSs in the bandstructure (II) with 
A = 0.05eV in the k z = n plane, (b)-(c) The obtained gap 
functions in bandstructure (II) in the k z = 7r-plane for (b) 
g = 0.127 eV and U = 0 (s++-wave state), and (c) U = 
0.434 eV and g = 0 (s±-wave state). The dots represent the 
experimental date given by ARPES study (in unit meV) [23]. 


LDA is much smaller than the reports by ARPES stud¬ 
ies. In order to magnify its size, we introduce the ad¬ 
ditional dependent potential T-Z line: Hkin = H° + 
zA/erfc e i\ (cos k x cos ky + l)cl{k)c l<7 {k), where e ; is the 
energy shift of the orbital l in the fc 2 -axis. We put 
e xy = 0.15 eV, and the others are 0. The obtained FSs 
in the k z = 0 plane are shown in Fig. 8 (a). Near the T 
point, the d^-orbital h-FS is largest of all h-FSs. Here, 
we neglect the SOI since the SOI-induced modification of 
the FSs is very small in Bal22 systems. 

Next, we solve the Eliashberg equation by introduc- 








































































10 



FIG. 8: (Color online) (a) The FSs of the three dimensional 
model for BaFe 2 As 2 in the k z = 0 plane. The green, red, blue, 
and black lines correspond to xz, yz, xy, and z 2 orbitals, 
respectively, (b) The obtained a c -a s phase diagram. The 
gap structure at each point a ~ c is shown in the figure. 
The horizontal axis is 9 = tan(k y /k x ), where (k x ,k y ) is the 
momentum on the FS with the origin at the center of each 
pocket, (c) ni m p-a s phase diagram in the superconducting 
state. The hole-s± state is realized for wide region in the 
phase diagram. 


ing dilute concentrations of impurities (1=1 eV and 
tiimp = 0.04) by applying the T-matrix approximation 
[18] to make the obtained gap structures smoother. Fig¬ 
ure 8 (b) shows the obtained a c -a s diagram of the su¬ 
perconducting state. The gap functions on the FSs (in 
the k z = 0 plane) at each point a ~ c are shown in the 


figure. The obtained phase diagram is qualitatively sim¬ 
ilar to that for LiFeAs shown in Fig. 5. Note that the 
gap structure in the region ”A!) y ~ 0 (A^, ~ 0)” is 

similar to the s++ (s±)-wave state. 

In the hole-s±-wave state at point b, only the gap func¬ 
tion on the outer h-FS (d^-orbital) is sign-reversed. We 
stress that the gap on the inner e-FS takes the minimum 
at 9 = 7 t /2 in the k z = 0-plane, which is consistent with 
ARPES measurement for BaFe 2 (As,P )2 [20]. This hole- 
s±-wave state is changed to the s++-wave state (point c) 
by increasing the impurities. During this change, large 
impurity-induced DOS appears at the Fermi level due 
to the outer h-FS, which may correspond to the exper¬ 
imental disorder-induced gap structure in BaFe 2 (As,P )2 
reported in Ref. [54]. 

Figure 8 (c) shows the ni mp -a s phase diagram in the 
superconducting state, where ni mp is concentration of the 
non-magnetic impurity (I = leV). The hole-s± state 
is realized for wide region in the phase diagram, re¬ 
flecting the fact that the hole-s± state is realized by 
the “cooperation” between orbital- and spin-fluctuations. 
That is, the “intra-orbital repulsive interaction by spin- 
fluctuations” mainly works between outer h-FS and e- 
FSs, whereas the “inter-orbital attractive interaction by 
orbital-fluctuations” mainly works between inner and 
middle h-FSs (d xz / yz - orbitals) and e-FSs. In the SC- 
VCe theory, both s ++ and hole-s± states are realized 
in the case of a c ~ a s even for rii mp = 0 , since the 
electron-orbiton coupling constant (~ U') is enlarged by 
the vertex corrections for the gap equation (A-VC) [53]. 
A crossover from hole-s± state to s++ state may be real¬ 
ized in BaFe 2 (As,P )2 by electron irradiation in Ref. [54], 


Appendix B: Spin susceptibility in the presence of 
the SOI 


In the main text, we studied the superconducting state 
in LiFeAs in the presence of the SOI. However, we ne¬ 
glected the SOI in the spin and charge susceptibilities in 
the pairing interactions V s (q) and V s (q), respectively, 
given in Eq. (16). This approximation is justified if the 
relation Xx(q) ~ Xy(q) ~ Xz(q) is satisfied, although 
this relation is not exact in the presence of the SOI. In 
order to verify this relation, here we calculate the spin 
susceptibilities in the presence of the SOI using the RPA. 

In the presence of the SOI, the spin is not conserved. 
Therefore, the Green function is given by the follow¬ 
ing 20 x 20 matrix form, G(k ) = (ie n + /i — h(k))- 1 , 
where h(k ) is given in Eq. (9). Then, the irreducible 
susceptibility is X^i'a'-^m'p'^) = G/ CT>mp (g + 

k)G m 'p' J'a' (^) ■ 

In the RPA, the susceptibility in the presence of the 
SOI is given by the following 20 2 x 20 2 matrix form: 

x(g) = (l-X°(q)f)- 1 X°(q), (Bl) 


where P is the Coulomb interaction in the matrix form: 




er,p' 


' o' ,p* 
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(a)/fe r .M_-=0,R)„=0) 


(b) / x (q x ,q y ,q^0,ai=0) 




a, = 0.98 
1=0 


a s = 
1 = 


0.98 

0.05eV 


180 




(d)lSfe,? r ?z=0,ra„=0) 



a, = 0.98 
1 = 0.05eV 



180 

120 


a, = 0.98 
1 = 0.05eV 


FIG. 9: (Color online) Spin susceptibilities in the case of as = 
0.98. (a) X s (q) f° r i n the absence of the SOI. (b) Xx(q ), (c) 
Xy(q ) and (d) X S z(q) in the presence of the SOI (A = 0.05 eV). 

(a) X(q„q y ,q=0,a>=0) (b) / x (q x ,q r q : = 0,aj„=0) 




FIG. 10: (Color online) Spin susceptibilities in the case of 
as = 0.90. (a) x "(<?) for in the absence of the SOI. (b) xJ(9)i 
(c) Xy{q) all d (d) Xz(q) in the presence of the SOI (A = 0.05 
eV). 


Then, the spin susceptibility in the /r-direction is 

Xfj,{q) = ^ ^ ' Xlcr,la';mp,m.p'(q) a 'a' ,cr a p,p' (B2) 

l,m,cr,a' ,p,p' 

where /r = x, y 1 z. 

In Fig. 9, we show the numerical results in the vicinity 
of the magnetic critical point, at which the spin Stoner 
factor is as = 0.98. Figure 9 shows the obtained spin 
susceptibilities for A = 0 (U = 0.439 eV) in (a) and 
for A = 0.05 eV (U = 0.446eV) in (b)-(d), respectively. 
In the case of as = 0.98, the effect of the SOI on the 
spin susceptibility is not small. In cotrast, we verified 
that the effect of the SOI on the charge susceptibility 
= Y,a,pXloM\mp,mp(<l) is very small even 
for ac = 0.98. 

In Fig. 10, we show the numerical results when the 
spin fluctuations are moderate (as = 0.90), which is 
consistent with the weak spin fluctuations observed in 
LiFeAs. Figure 9 shows the obtained spin susceptibilities 
for A = 0 (U = 0.403 eV) in (a) and for A = 0.05 eV 
(U = 0.409 eV) in (b)-(d), respectively. In the case of 
a s = 0.90, the relation y*(q) ~ x£(<?) ~ xKq) is ap¬ 
proximately satisfied, so the rotational symmetry of the 
spin susceptibility is almost satisfied. Thus, the effect 
of the SOI on the pairing interaction is expected to be 
negligible. 

Experimentally, LiFeAs is not close to the magnetic 
critical point because the observed spin fluctuations are 
small, so the relation as ~ 0.9 is expected in LiFeAs. We 
have verified that the gap structure in the s-t-state shown 
in Fig. 3 is essentially unchanged even for as = 0.90, 
similarly to our previous study in Ref. [32]. Therefore, 
the effect of the SOI on the pairing interaction V s,c (q) in 
LiFeAs is negligible, so the obtained results in the main 
text are essentially unchanged even if the SOI is included 
in the pairing interaction. 
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